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ABSTRACT 


Atmospheric density fluctuations induce phase 
perturbations that degrade the spatial coherence of a laser 
beam propagating through the atmosphere. These 
degradations spread the laser beam and alter the centroid 
and intensity profile stochastically. Turbulent conditions 
are found at virtually all levels of the atmosphere. A 
substantial fraction of the optical turbulence along a 
vertical path arises from the heat flux between the 
atmosphere and the Earth’s Surtace. This type of 
turbulence is typically within the first 100 to 200 meters 
above the surface. 

During this thesis research, a high frequency acoustic 
sounder was developed to analyze this turbulent layer. The 
primary focus was the development of the command and 
control software required to coordinate the data collection 
and reduction. The system was used at two sites and should 
prove useful in quantifying the eftects Of optical 
turbulence within the surtace boundary layer on laser and 


optical system performance. 
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I. INTRODUCTION 


A coherent laser beam propagating through the 
atmosphere iS very susceptible to numerous turbulence 
inflicted degradations. As electromagnetic waves transit 
the turbulent regions, atmospheric irregularities randomize 
the amplitude and phase of the wave. In order to quantify 
the altitude dependence of the atmospheric turbulence, a 
high resolution turbulence profiler is needed. Acoustic 
echosounders are frequently used to detect and measure 
atmospheric density and velocity irregularities resulting 
Erom air currents, temperature inversions, humidity 
variations, mechanical turbulence ano s OLnerm Censluey 
fluctuations. 

Presently two atmospheric optical parameters, the 
spatial coherence length (ro) and the isoplanatic angle 
(80) are measures of the perturbation of an electromagnetic 
wave propagating through the atmosphere, and are accurately 
measured by optical systems developed by Walters [{Refs. 1 
amd 2 |e Although these systems measure a path integral of 
atmospheric turbulence with high accuracy, a major drawback 
of each system is that no provision is available for 
determining the height of the atmospheric disruptions. it 
these disturbances are found to exist very near the 


surface, it may be possible to negate their effect by 


elevating the optical systems or controlling the generation 
of turbulence: 

This thesis deals with the design, construction and 
implementation of a high frequency, acoustic echosounder 
which will accurately analyze the atmospheric density 
fivetuationus © wien approximately 200 meters of the 
surface. As this project is a product of the research and 
efforts of two students, the work was appropriately 
divided. My particular task was to devise a computer 
program which would control all aspects of the echosounder 
operation. Areas’ of particular interest were the 
integration of hardware and software within the data 
acquistion system, developing ame 1a oS algorithms, 
controlling the system timing, and setting the input 
parameters which determine the range and sensitivity of the 
device. My colleage, LT. Weingartner [{Ref. 3], dedicated 
his efforts toward the actual design and hardware 
development of the echosounder. 

Acoustic echosounders have been developed and in use 
for many years and have proven to be valuable probes for 
analyzing the structure and dynamics of the lower 
atmosphere [Refs. 4 through 7]. Devices similar to ours 
have been used to obtain profiles of the atmospheric 
density and temperature fluctuations [Refs. 8 and 9]. Our 
device used a high speed HP217 computer to control and 


monitor the echosounder which provided real time 


tnitormatvion. Additionally, we had the ability to store and 


reproduce the atmospheric profile plots at will. This 
information was crucial in analyzing the atmospheric 
turbulence degradations on laser and electro-optical 


system’s performance. 


This document shall address the theoretical background 


of acoustic sounder operation in Chapter II. Chapter III 
shall present a synopsis of the system design and 
associated software. A summary of actual data collected at 


two sites is analyzed in Chapter IV and Chapter V discusses 
the conclusions and recommendations found as a result of 


our research. 


Il. BACKGROUND 


Acoustic echosounders probe the atmosphere by 
transmitting a pulse of acoustic power which is 
subsequently scattered back from the atmosphere by 
temperature and velocity inhomogeneities. The echosounder 
(echosonde) equation, often referred to as the radar 
equation in meteorology, is used NS determine the 
backscattered acoustic power. This equation is summarized 
by Neff in Reference 10 and is based upon the work of 


Tatarski {Ref. 11] and Little [Ref. 5}. 


Pr = ErR{ PrEr }f{ exp(-2ak) Ji coCRee jet seu |) “Ghee 


where 

PR is the electrical power returned from a range R, 

Pr is the electrical power transmitted at frequency f, 

KR is the efficiency of conversion from acoustic power 
to electrical power by the transducer, ° 

Et is the efficiency of conversion from electrical 
power to acoustic power by the transducer, 

exp(-2ak) is the round trip power loss due to 
attenuation where ais the average attenuation (meters-! ) 
to the scattering volume at the range R (meters), 


oo0(R,f) is the scattering cross section per unit volume 
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at a distance R and frequency f, 

ec is the local speed of sound (meters/second), 

t is the pulse length (seconds), 

A is the aperture area of the antenna (meters® ), 

R is the range (meters), and 

G is the effective-aperture factor of the antenna. 

Empirically measuring or calculating the values for all 
other terms, one can use this equation to determine 
co(R,f), the scattering cross section per unit volume; that 
is, the fraction of incident power backscattered per unit 
distance into aunit solid angle at a frequency f. Based 


upon experimental results, Tatarski [Ref. 11] expresses the 


acoustic backscatter cross section per unit volume, 
Some: js in the equation, 
Cre 
oo0(R,f) = 0.0039 ki/3 ceenrnen- : 
EO 
where 
Kee 2s A is the incident acoustic wavenumber at 


wavelength i, 
To is the local mean temperature in degrees Kelvin, and 
CtT* is the temperature structure parameter. 

Combining this equation with the echosonde equation, one 


obtains a volume-averaged measure of CT®. 


Ld 


i i Zs 1 PR 
Ci? 3 === == 2-222 To* k- 1/3 333 =  s= >= == he ariel! 
0.0039 wren Hr ec AG Est 
Hall and Wescott [{Ref. 12} calculated a beam-shape 
compensation factor of 0.40 for a piston source antenna 
with a uniformly illuminated square aperture. This value 
is the same as the effective-aperture factor, G, and can be 
substituted into the above equation. Approximating the 
aperture area of the antenna to be equivalent to 25 times 
the aperture area of a single speaker having a diameter of 
7.620 centimeters, we get a value of 0.1140 square meters. 
Combining these values with the numerical constants in the 


above equation, we can simplify the equation for CtT?. 


it 1 PR 
Cr* = 11245 ------ To? k-1/3 eee- ---- Re exp( 2aR) 
ER Et CT Ey 
The acousto-electrical efficiency factors for our 


echosounder were measured with a calibrated microphone and 
an anechoic chamber. The calculated values for ER and ET 
are each 0.5. Substituting these values into the above 


equation, we can further simplify the equation for CT*. 


To == GER 
Ct? = 44980 k-1/3 e2-2-2- ---- R? exp(2ak) 
CG Ea 


iz 


The above equation was incorporated into the computer 
program that provides the system control, data acquistion 
and data processing techniques [Appendix A]. The- reduced 
data was then used to measure the temperature structure 
parameter as a function of time of day and altitude for 


various sites. 
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SYSTEM DESIGN AND EQUIPMENT DEVELOPMENT 


used in this system was standard 


and is illustrated in Figure 1. The 


mpvyvcnoroz2zm 


CAAT 
5 X5 Speaker 
Array 


Echosounder Equipment Arrangement 
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pieces which were 


‘pecially designed by LT Weingartner 
Peers are the acoltliscic array and the enclosure for this 
array. The design considerations for the speakers’ and the 


Array format are given in Appendix B and the desi the 


ay 
= 
a 


enclosure to house this array is outlined in Appendix C. 


Aside from these two pieces ot equipment, there are 
primarily Six other components which comprises this 
echosounder. Below is a brist description of each ot these 


additional components. 


Ll. HP200 Series Computer 





The Hewlett-Packard (HP) 200 Series Compnter 
includes a 20 megabyte hard disk drive, a tloppy disk drive 
and an associated printer and monitor. ine HPAZ0O Series 
Computer was the central control component for the entire 
echosounder arrangement. The computer used was an HPY1; 


Peowewommea an Basic 3.0 and equipped with an Infotek BCUUs 


posi  toMmp. Ler and an ioe teeee Oo Blaating “Faint 
Accelerator to enhance the speed of program execution. Thre 


program detined all parameters for the HP3314A Function 
Geeeeator as well as execute the trigger command which 


produced the echosounder transmitted acoustic pulse. The 


computer also received the data from the acoustic array Viz 


the pre-amplifier, bandpass filter and analog to digital 
eonverter. The computer then conducted the data reduction 
bev ines Lo produce ama display ~a high resolution 


al! 


BMOSpueric pretile. 


2. HP3314A Function Generator 
The Hewlett-Packard (HP) 3314A Function Generator 
was a multimode function generator capable of providing 
Sine, square and triangular wave functions as well as any 
desired wavetorm ranging in frequency from 0.001 Hertz to 
19.999 Megahertz. The HP3314A Function Generator was used 
to supply a pulse of an integer number (usually 100) of 
Sinusoidal cycles of constant amplitude to the QSC Model 
1700 Audio Amplifier. <A constant frequency setting of 5000 
Hertz was used for all data runs. 
3. @QSC Model 1700 Audio Amplifier 
The QSC Model 1700 Audio Amplifier is a high power 
amplifier which can supply 350 watts over an 8 ohm load. 
This amplifier was used to boost the output voltage of the 
function generator By a factor of ZOMinenme ovo !ts tome 
Vougsce The power supplied to each speaker in the array 
during operation was then 37.5 watts. 
4. Pre-amplifier 
The pre-amplifier was designed and constructed by 
Walters [{Ref. 13} to supply a gain of 10005 to the returned 
Signal and isolate the data acquistion components from the 
transmitted pulse. The pre-amplifier may be thought of as 
a safety and switching mechanism for the system. An LT 
1037 and an LT 1007 Operational Amplifier were selected for 
use in the pre-amplifier based upon their low noise 


properties which were evaluated by physically incorporating 
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them into the device and measuring their noise and gain 
characteristics with an HP3561 Spectrum Analyzer. The 
diodes used are IN 4000 Series Rectifiers capable of 
carrying one ampere. These rectifiers serve to isolate the 
power amplifier from the amplifier input stage and to limit 
the voltage applied to the operational amplifiers. The 


electrical diagram of this device is outlined Figure 2. 


5602 TO ARRAY 
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Fig. 2. Schematic of Pre-amplifier 


ie 


5. Rockland Wavetek Model 852 Filter 
The Rockland Wavetek Model 852 Filter operated as a 
48 db per octave bandpass filter to suppress the broadband 
noise of the system. High and low bandpass settings of 
5500 Hertz and 4500 Hertz respectively were used. The 
filter response at these settings is illustrated in Figure 


on 
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Fig. 3. Bandpass Filter Response 


The three db (half power) bandwidth is approximately 1 KHz 
wide. Although it would be desirable to reduce the 
bandpass to between 50 and 100 Hertz without impairing 


Signal response, the Rockland Wavetek Model 852 Filter was 
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the best tilter available. It performed adequately in 
spite of the large 1KH2 bandwidth. 
6. Intotek AD200 Analog to Digital Converter 
Mii ntouemeAuUZzOUml2 BiG Analos to Digital Converter 
was used to digitize the signal voltage for the computer 


system at a 12.5 to 20 KHz sample rate. 


B. SOFTWARE 

The software and HP200 Series Computer system were 
responsible for controlling and monitoring every 
operational phase of the hardware components. The program 
which accomplished this task was entitled "ACRDR" and is 
listed in Appendix A. This program was modeled after a 
program written by Walters [Ref. 13] but each program 
performs a distinctly different computational task. nas 
program was written in HP Basic 3.0 and was compiled by an 
Infotek BC203 Basic Compiler to enhance the speed of 
execution. 

The program "“ACRDR" is easily broken down into a number 
of blocks and subroutines which performed specific 
operations. These sections are outlined in a flowchart 
(Fig. 4). Each block is straightforward in its purpose, 
and the program is designed to be as helpful to the user as 
possible. As a prologue to the actual program code, there 
is a listing of all the program variables with a short 
description of their use. Such a listing familiarizes a 
user with the computations to be made and also provides 


i” 


quick references for any ofuturemenodificattonsmet ome 


implemented. 
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Fig. 4. Flowchart of Fundamental Program Operations 


The initzwalization procedures) which set all the 
parameters necessary for data collection follow the 
variable definition section. These parameters control such 


things as the contrast between background and return signal 


and the setting of the computer’s internal clock. All 
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options of the HP3314A Function Generator are decided by 
the user at this’ phase. Another initialization involves 
preparing the Infotek AD2ZO0OO Analog to Digital Converter for 
operation. Finally, internal arrays are dimensioned, 
integer variables are defined and the computer function 
keys are redefined to suit the “ACRDR” program. 

Now that the system is ready to operate, certain 
quantities must be determined based upon atmospheric 
Cahicdtirsems. atwathe  stime of date collection. From the 
ambient temperature, a speed of sound in air is determined 


by the equation in Kinsler, et al, [{Reft. 14: p. 106] 


te Cont fe Cl72ta) yi7e , 


where 

C is the local speed of sound, 

Co iS the speed of sound in air at 0” Celsius (331.6 
m/s), and 

T is the local temperature in degrees Celsius. 
Using this speed with a user decided sampling frequency and 
number of data points, the maximum range of the traces 15s 


determined by the relation, 


RANGE = C/2 ¢ 1/(SAMPLE FREQUENCY) ¢- N , 


where 
RANGE is the range of the echosounder in meters, 


SAMPLE FREQUENCY is the user desired rate of data 


Zu 


acquisition, and 

N is the number of data points per trace. 

Typically, the values of many of the parameters are not 
varied. The value of N above is almost always taken to be 
the allowed maximum of 16301 data points and the speed of 
sound can be safely estimated to be 340 m/s at room 
temperature. By using a typical sample frequency of 20000 
Hz, a range of approximately 135 meters is obtained. 

Another parameter determined is the atmospheric 
attenuation coefficient. This calculation is done ina 
subroutine obtained from Reference 15 and is necessary in 
this program for the calculation oF Cre2. Other 
computations made prior to program execution include the 
pulse length, wavenumber, and the D-C offset of the 
equipment. 

At this point, the computer is finally ready for data 
collection. The screen setup displays a distance versus 
time plot and the internal clock of the computer is 
synchronized with the pulse from the HP3314A Function 
Generator. Following transmission of the pulse, the return 
Signal is received, digitized, stored in an array and data 
reduction commences. A block averaging technique is used 
in which a block of data points the size of the number of 
emitted Benes ae summed and averaged. Returns are plotted 
as darkened areas on the echosounder traces with the 


intensity of the darkened area linearly related to the 


ee 


Magnitude ot the return signal. This same procedure is 
used in the computation of CT2 with some very important 
differences. Aside from simply correcting tor nr 
spherical divergence, CTt2 is corrected for electronic gain, 
the ratio of power returned to power transmitted, the 
efficiency of the speakers as transmitters and receivers, 
the area of the speaker array, the atmospheric attenuation, 
the pulse length, the temperature, the scattering cross 
section per unit volume at a specific range and frequency 
and finally the effective aperture factor of the antenna. 
Additionally, Cr2(Range) is averaged for a particular 
altitude over 15 minute intervals. Finally, after each 
pulse is reduced, a corresponding mean square noise level 
is determined. This is done by averaging each block 
average at maximum range until at least ten values have 
been used in the average. An upper limit of five over the 
average noise figure (nies corresponds to voltage 
fluctuations on the order of 10-7 volts) is set on the 
routine to avoid averaging any strong return signals or 
anomalies such as passing aircraft. Atter the ten values 
are averaged, every subsequent pulse is averaged into all 
the preceding noise levels and removed from each subsequent 
return signal. 

Ultimately, after each 15 minute interval, or at the 
users request, the display terminal image is printed. PSG 


this time the Ctr2 computations are conducted and plotted on 


ae 


the screen and printer as range versus logio of Cr2. Upon 
completion of the printing of these plots, the system 
begins the data processing for the next 15 minute interval. 

There are certain options built into the program to 
allow a user to change various aspects of operation. The 
function keys allow the user to change the local 
temperature, sample frequency or intensity factor during 
program execution. Additionally, the user can quit or 
restart the program, print the partial trace on the screen 
or elect to save a trace on a floppy disc. The save 
routine 1s invoked for the subsequent 15 minute interval 
after the appropriate function key is depressed. Saving a 
future trace may seem awkward, especially if a user would 
like to keep an interesting trace which is presently on the 
terminal. This problem cannot be readily solved unless 
each trace is recorded to disc without user intervention. 
At present, this is not done because only 8 traces (Z hours 


of data) can be written to a floppy disc before it is full. 


IV. DATA ANALYSIS 


A.  ECHOSOUNDER PERFORMANCE 

An analysis of the acoustic echosounder output was 
conducted to determine the validity of previously 
determined echosounder parameters. The e-1 decay time 
constant for applied voltage was calculated in Appendix B 
to be approximately 900 us or for convenience, 1 ms. 
Reviewing typical echosounder traces and Ct2 plots 
[Appendix D}], it was evident that the recovery time was 
typically found to be on the order of 33 ms or roughly six 
meters past the end of the transmitted pulse length of 3.4 
meters. The recovery time iS consistent with the time 
required for the 30 volts on the drivers to decay to the 
microvolt level. A more detailed analysis of the hardware 


is found in Reference 3. 


B. SITE EVALUATION 

Rchosounder data was collected at two different 
locations. The primary data collection site was the upper 
roof of Spanagel Hall at the Naval Postgraduate school, 
Monterey, California. This site was chosen simply for 
convenience. Data gathered at this location is believed to 
represent the California coast during the spring near sea 
level. The second site chosen was in the vicinity of the 


2 Ae telescope at Lick Observatory, San Jose, 


Za 


California. This site is lo@ated atepeit.) dani ene ee 
altitude of approximately 5700 feet and nearly 20 miles 
inland from the coast. 

These two data collection sites represent areas of 
differing atmospheric air pressures, water vapor pressures, 
local temperature ranges, and local wind velocity ranges. 
These characteristics all play important roles in effecting 
the local atmospheric turbulent conditions and thereby the 
atmospheric structure parameter, CtT2. 

In addition to collecting echosounder data at Lick 
Observatory, simultaneous measurements of the isoplanatic 
angle (80) and spatial coherence length (ro) were made with 
systems developed by Walters [Refs. ierana s 2 | A basic 
knowledge of these two systems is necessary to understand 
the correlation procedures’ made. The isoplanatic angle 
(60) is primarily an upper atmospheric measurement which 
indicates atmospheric disruptions at a range of 2 to 15 
kilometers. The spatial coherence length (ro) is a measure 
of the effects of the entire atmospheric blanket on 
coherent light transmission. A close comparison of all 
three data sets should give us an accurate description of 
both the lower and upper troposphere as well as the 
stratosphere above. 

Based upon isoplanatic angle (80) and spatial coherence 
length (ro) measurements at Mt. Wilson in California, a 


strong correlation between the two measurements occurs if 
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the low altitude boundary layer COnvr Duct on is 
Sime eetedmy smal 1. This strong correlation helps to 
reinforce the overall description of the atmosphere at the 
time of data collection. Figures 5 and 6 EEO aie 
illustrate the atmospheric measurements made at Mt. Wilson 


on 2 April 1987 by Walters. 
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Fig. 5. Isoplanatic Angle Measurements, Mt. Wilson 
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Fig. 6. Spatial Coherence Length Measurements, Mt. Wilson 


The strong correlation between the isoplanatic angle 
(60) and the spatial coherence length (ro) is especially 
evident between the hours of 0700 and 1300 universal time. 
The close tracking of these two measurements during this 
time interval indicate that the upper atmospheric 
conditions, as measured by 80, are dominating the entire 
atmospheric profile as measured by ro. Unfortunately, the 
existence or non-existence of any turbulent surface effects 
cannot be ascertained by the employment of the above two 


systems alone. However, employment of these two systems 
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together with the acoustic echosounder should enable us to 
produce a complete atmospheric profile with strong 
correlation between all three atmospheric measurements. 

On 9 and 10 April 1987, all three systems were operated 
at Lick Observatory. Again, a good correlation between the 
isoplanatic angle (80) and the spatial coherence length 
(ro) measurements was noted. In addition, a strong 
correlation between the spatial coherence length (ro) and 
echosounder measurements was present. A comparison of the 
atmospheric data in Figures 7 and 8 shows a good 
correlation between the two parameters especially during 
the 0930 to 1130 time interval on 9 April 1987. However, 
during subsequent hours the isoplanatic angle (80) 
measurements remain high (=12 urad) indicating relatively 
calm turbulent conditions in the upper atmosphere while the 
spatial coherence length (ro) values drop sharply after 
12:00 Universal Time indicating dominant and increasing 
lower atmospheric turbulence. This increase in the low 
level turbulence should be evident in the echosounder data 
commencing around 1200 universal time (0400 local standard 
time) on 9 April 1987. A comparison of the echosounder 
data in Figures 9 through 11 illustrates this increase in 


the local surface turbulence. 
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Fig. 8. Spatial Coherence Length Measurements 
Lick Observatory 
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Fig. 11. Echosounder Trace, Lick Observatory 


The strong correlation between the echosounder data around 
4:00 to 5:00 Standard Time and the ro measurements around 
12:00 Universal Time combined with the lack of correlation 
between the 980 and ro measurements indicate that the lower 
atmospheric and surface turbulence are dominating the 
atmospheric profile during this time period. 

Data collected on 10 April 1987, again illustrate the 
strong correlation between the three atmospheric 
measurements’ made. A comparison of the _ 60 and ro 
measurements during the time interval of 0700 and 1300 
Universal Time indicate steady turbulent conditions in the 


upper atmosphere and greatly varying turbulent conditions 


rz 


at lower atmospheric levels. This is evident in Figures 12 
and 13 by the consistent values of 80 during the time 
period compared with the steady increase and eventual 
decline of the ro values during the same time interval. 
The trace variations in Figure 13 during the hours of 0800 
and 1200 are indicative of a period of decreasing lower 
atmospheric or surface turbulence followed by the onset of 
an increasingly turbulent period around 1200. This 
turbulent trend is strongly supported by the echosounder 


data in Figures 14 through 21. 
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Fig. 12. Isoplanatic Angle Measurements, Lick Observatory 


Sica! 





LSS Zee leaiee. 


1G Aerial 


36a 


10 le Fa 


Chiesa) 


8 
Time 


13. Spatial Coherence Length Measurements, 


Fie. 


Lick Observatory 


LICK GbSeRyarorRy 


L267 


10 Apr 


Tso 





Sd3LAW NI ADNUY 


Tk Et OC Gee S aD 


Lick Observatory 


Echosounder Trace, 


14. 


Fig. 


34 


So ae 
Y .. 
oe 
reas 

ie 
a 
tJ ...: 
UW) , 
ny . 
O oO 
md 

O 

| 

—) 

f 

CO 

Oo) 

< 

Q 

oa 

Q 


135 


4m: 


iS 


1:18 


PEE Oear  (s1D) 


BS 


| 





OG 


SY3LaAW NI SONYA 


Lick Observatory 


15. Echosounder Trace, 


Fig. 


PICK Greer VATORY 


1@ Apr 198? 


eae 





oe eset Ay 
: i wa tg bat 
o fl ae 6 0 Sse 
. e . od 
oF ee 
: weg 
a “ie e syd 
ae e: : wits ee 
oe! ® . a “eo : er ad 
ABO BP 408 ei 
*% i a¥ 
. “Rat el “Ty 
e 3 eG! fs ¥ = 
w. Loo ancetsag A et a 
Mom peeesearere ears “SoM Ses) | eae t 
ee ee: Meserti ee 2 PURE] 
‘we atte ee ree am 
Mea aan mes EM 
. ate . uN 
C c oe ea “> 
ee ‘~ ¢ 8 e = a > gt 
- = ee oe umenstee® ne : 
e « : 6 gts, hres ai A : Us 
* . . Hr ] 
te He ee : ee nh at ae 
Ses : La ey arte dij 
. oJ * a> 
« a oe nant ma vrai 1 
oe t+ ee 8 
we . ° Dw Pye 
e . 5 a, oe : 5 
Law  a TE HA 
eee aes gee a ; ns 
peeees one be Th A) 
o Me Habe SS a E 
. o” oe STMAR USS e — 
oe : rie ga ee vs; f 
ia © RE eS Re A rg es DR 
cee Fe OUR ee cee ee set! hie 
on © one of. CE aeons : . re a “ 
gifs. 4 rf "y 
‘Seay i 
: Hi! aval lobe, SO 
- . 2 TEP Sa 
oxen ves PNT 
. ws Me 
e es a et EVE 
ene oe e¢ @ 5 -. Re yee ws edgy i 
er ceneasrcas . rye 
‘ = SE ah ae 
yea oR : ute y 


is 


SYILAW NI 3DNHY 


Pe eeOe Rin Cs 1) 


Lick Observatory 


Echosounder Trace, 


ioe 


ae 


oe 


RANGE IN METERS 


RANGE IN METERS 


30 


Tee Reeee sie LICK (Gees FiIORy 






=) e 2 we 
° . e ° oe 8 . 
ne e re ‘ . 5 ° . j ° ° s ORS . - 
ie . - ° @:¢e es eee e #¢ w 3 
. . e ee . ee e ee . e ° ‘ Let 
. . ee e a ° 
ne Sag oe ese a 
a * oon e * ee . r 
. : -_«” . 
. = vo g ° Tie wer 
e s 







. 
ome wunene 


2:22 2:05 2:18 15 


TIMe ECGhe Coie 


Fig. 17. Echosounder Trace, Lick Observatory 
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Fig. 21. Echosounder Trace, Lick Observatory 


The large 200mm coherence lengths around 11:00 Universal 
Time are consistent with the low turbulence evident in the 
echosounder profiles around 3:00 Standard Time followed by 
a pre-dawn increase in the surface turbulence. 

Only acoustic sounder atmospheric measurements were 
made on the upper roof of Spanagel Hall at the Naval 
Postgraduate School. Data runs on 26 and 27 April 1987, 
included both echosounder data and the associated Cr? 
plots: A representative sample of the data collected 
during this period is included in Appendix OD. Of 
particular interest is the data in Figure 22, which shows 
the maritime boundary inversion layer at about 100m and 


being perturbed by convective plumes at lower altitudes. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A high frequency (5KHz) phased array echosounder 
constructed to measure low level turbulence appears to work 
within the 10 to 135 meter altitude range. When using the 
device, detailed profiles of the short range atmospheric 
density fluctuations were obtained. The profiles were 
found to correlate very well with the measurements of the 
isoplanatic angle (80) and the spatial coherence length 
(ro) during periods of simultaneous operation. This short 
range echosounder, when used in conjunction with the other 
atmospheric measuring devices, is an invaluable tool. The 
graphical output provides us with a more complete 
description of the atmosphere that can be used to calculate 
Chae the atmospheric index of refraction structure 
parameter. 

Further software development such as the incorporation 
of Fast Fourier Transform routines into the echosounder 
program will determine the radial velocity profile of the 
return signal. such an improvement will determine the 
velocity of the probed air masses and plot these 
echosounder traces as a function of color intensity. 

Other areas of further reasearch include the design and 
testing of different array patterns. One such pattern, the 


hexagonal array, is already undergoing tests. As the 


40) 


speaker array evolves, it is also evident that the 
enclosure must follow suit to accomodate the new array 
design and associated beam patterns. Further improvements 
in software routines are also inevitable. It is always 
desirable to store on a dise all data collected at a site. 
Presently, the floppy disc capacity is inadequate for the 
storage of more than two hours of data. Furthermore, if 
FFT and Doppler routines are to be added, the computational 
speed of the computer will inhibit the pulse repetition 
rate. A special data acquisition technique called “Direct 
Memory Access" may then have to be introduced into the HP 
computer to allow Simultaneous data collection and 
processing. 

Haemakl y , hes DrOauce. may be used in various 
applications not already addressed. One such use may be to 
measure windshear at airports. Other researchers have 


expressed interest in using the echosounder to measure 


eye ce atmospheric conditions during meteorological 
surveys. Assuredly, as the device evolves into a more 
compact and highly efficient instrument, its range of 


application will continue to expand. 
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APPENDIX A 


ACOUSTIC ECHOSOUNDER PROGRAM 


KEEKKEHREKEEEEEEKEEEEEEEEEEREEEEE EERE EEE REESE EEE 


i ACRDR 

* COMPUTER SOFTWARE WRITTEN BT LT. M. WROBLEWSKI FOR AN 

* ECHOSGUNDER BUILT FOR Aetiasteks CF SClENCe (hes 

BY LT. M. WROBLEWSKI ANO LT. F. WEINGARTNER 

* ADVISOR ERO =e eweren > 

% JUN. 1987 
EEEEKEEKEEKEKEEEEKEEKEEEKEEEEEKEEEERKEEEEREKE EKER ESE 


The computer program receives information from an acoustic 
array through an A-D converter. This information is the 
returned signal of an acoustic pulse as it passes through 


the atmosphere. The data is then used to display the return 


intensity with distance as a function of time. This is a 
short range device (from 158m to 20@m). 


ELST Geaverlogiees 
I& kK & J - counters for loops 


Rec_num & Nrec - real and integer representation of the 
record number for storage 


Plotnum - counter used to insert form feeds between plots 

Disc_adcdressS - storage location of data file 

FileS - filé mame UsecG tao store data 

NS &@ Num - string and réai representation ef the number of 
cyclés and the numter of points used in the 
computation of the Glock average 

Point(*) -. réal number representing the noise and range 


corrected average over given number of points 


Dat(*) - array of a-d converter output after sampling 


42 


* 
* 
¥ 
* 
# 
¥ 
Lg 


Cd Cd Cd G3 GD CE 


O) oot) ©) Coo. c=] 
OF U1 t= Gra -- O13 CG 
G> C2 63 


ao @ oO 
co @M ~] 
ao Oa @& 


999 


Gaetemee- array Cfereeduced and averaged) data including 
offset, correction for noise, and range 
corrections 


Hesee — the integer hours 
Min - the number of minutes 
OCGnGeevOtemime —- Keeps track of the passage of each 


15 minute interval 


Nem! —- counter for conputation of block™averages 
Num2 - counter for the number of Block averages made 
M - the counter used in the plot label routines 
Timent —- time interval Between data samples in 
nanoseconds 
Kmax ~- total mumber of block averages computed 
Pressmoerceo — themrequencyminogt of the HP3Z314A in kilohertz 
ZoneS - the appropriate time cone the operator desires 
Slot le— ties heme Oh th appropriate "Siteof date collection 
Save _plt & Cntrl - on/off toggles used to determine whether 


a particular run will be saved to a disc 


Nolot - the number of increment along the entire 
horizontal axis 


Offset - the computed O-C offset for the system prior 
tO data acaguisition 


Noise - arunning total of the noise accumulated at the 
MAaXLMUM TENGE; UsSeT tO Tins ens ceaove the 
ave: auc horses 

Limit - the upper bounc on the noise figure used to insure 


that large returns are not included in the 
noise computations 


Ki) - a counter of the number of traces used to determine 
the noise average 
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Samave - a sample average used in the computation of the 
O=€ oftset ? 


Off_ave - a toggle used to inhibit the firing of the 
HP32Z14A during D-C offset computations 


CountS & Cnt - the string and real representation of 
the number of data points to be taken 


Samfreq - the user input sampling frequency desired 

Tminc - 1 over the sampling frequency; the time between 
sanples 

TimeS - the string required by the A-D converter to 


sample at the desired rate 


COkelvin - the speed of sound at @ degrees celcius in m/s 

Temp - temperature of the surroundings in degrees celcius 

Spd —- the speed of sound in air computed for the input 
temperature 

Rspd - the relative speed of sound that the echosounder 


sees which is half the computed speed 
Rdist - the distance traveled in one time increment 
Maxcng - the maximum range of the echosounder found 


by multiplying the distance per time increment 
by the number of time increments 


Exrng - the range rounded to the nearest value evenly 
divided by 15; used solely for plotting 
purposes 

Far - used to estimate the far field; this value keeps the 


range correction from being applied to data very near 
to the source 


Iivl - the user input intensity lever Givise sethiscme oluc 


seis the screen contrast in data return 


Scale_y - the ratio of éxrng to maxeng; this value is used 
to keep the plotted range of data accurate 


Npoint & Npoint! & Npoint2 - variables used to keep track of 
time passage between pulses 
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meme a — —a Ce ee ee eS 
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Cd a) 


Cease ier i109 — “Used to synchronize data collection with the 


clock 
(enol tee wet neewmternal Clock Of the comouter 
Ovd - divisor of the data average; used because computer 


multiplication is faster than division 


R - the range of a block of data samples 

Run_ave - the running sum of the block samples 

X - the value of the data point less the D-C offset 

Ns - the final average of the block of data points; also 


used in the computation of the noise figure 


CSgeee-metme moise cormPection applied to the data samples 
Time - the horizontal position of the trace on the plot 
Timedist - the horizontal width of the trace on the plot 
Endpt - the last point of the plist verticaily taking 
into account the Scaring factor scale_y 
Inc - The vertical increment along the plot 
Z - The final réduced data points which are output on the plot 
Cis - The vertical height of the trace on the plot 
Aas Fras ,NmS ,EnS Vos ,H2S5 - strings needed to set the HP3SI<AA 
AmpS - userv input amplitude for the HPS3SI4A 
Pls_lng - the pulse length of the burst 
Ct(+) - the atmospheric temperature structure parameter 
i. SemmumMmoee gt tiaces Seca im computing Ct. *) 
THM cnclvecticag! Dositlon on the Ct pict 
PGi > tme mMeGizomtal Gositicon om the Ct pleat 
KZ - the wave mumber to the I/3rd power 
A - The area of the receiver (array) 
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| G - The effective aperature factor 

| 

|! Er =~ efficiency of coversion of acoustical power 
| to electrical power on the recieve side 

! Et - the efficiency of conversion of electrical 
| power to acoustical power on the transmit side 
| 

! Pt - The computed transmitted power to the 
acoustic array 

| 

| Gain - the electronic gain of the equipment 

| 

| Zimp - the speaker array impedence 

| 

| 

OPTION BASE ! 


— — om — — oe 


initialize the arrays % 
Geclare all integer variables 


set dimensions 


DIM Disc_address3(20] ,FileS(30] Point( 320) ,Ct( 200) 
INTEGER I ,Hr ,Plotnum ,Print_key ,Num ,Numt ,Num2 MK 
INTEGER Rec_num ,Kmax 

INTEGER 02(300) ,Dat( 16301) BUFFER 


mM — eee oe 


R 


— oe rams sae — ° on 


initialization routines....sat time, sat HP3Z314A Function 


Generator 


tet :CaAlle Prégmimii. Preece. Ge 


CALE Init vadZoo 

CALL Set_time( Zoned) 
INPUTS ss ew NGNe seca. 
Lt=2 


SaVE_fB 
Cntirli=s@? 
jiu 27m} 


cose oe & 


Ss labels on the computer function 


OUTPUT KBOs “SCRATCHRie (om, LSGEEAR” THE skews 


CONTROL 2,231 
ON KEY 1S GRBEs 


"NEW TEMP" GOTO Speed 
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ON KEY 2 LABEL “INTENS. BocvOR., GOlrOrLVI 
ON KEY 4 LABEL “PRINT TRACE” GOTO Prt_dmp 
ON KEY S LABEL “SAVE NXT PLOT” GOTO Svpt 
ON KEY 6 LABEL “CHANGE SMPL FRQ" GOTO Sfrq 
ON KEY 7 LABEL “RESTART” GOTO Rstrt 
Cieiwerec LAGEL “QUIT GOTGsGui + 

Set constants 


Disc_addressS=": ,700,1 ,0" 


Maxrec 
Nolot= 


=(20 


920 


Plotnum=1 
Of Ffsat=2 


GINIT 


Noise= 
Limit= 


NS ome 


Kl=@ 


Q 
iQ 


Samave=Q 
A=.1142 


G=.4 


Er=.496 
Et=.496 
Gain=723126"%< 


Zimp=t2.12 


| Input and calculation of terms used by the A-D converter and 
the computation of range 
| 


PRINT 
PRINT 
PRINT 
PRINT 
INPUT 


1 Aan 


"Enter the number of data points desired (max 163901)" 


"Use increments of the number of cycles then add |” 
"ex: (150 cycles K 100 data points/cycle)+1=I1500i" 
"RECOMMEND USING MAX VALUE OF 162301" 

Counts 

BUt Power CoMmMButec 4s5s5uNEd 

plifier settings 


PtosetecoO a) )/(c*Zime) 


PRINT 


seen 


PRINT 
PRINT 
INPUT 


This frequency will" 
Examples are:" 
225 


oo 


"Input the sampling frequency desired. 
"determine the range of the sounder. 
i See REO me | 2500 wees anc RANGE 
E SAM. FREQ. = 20000....... RANGE 
"ENTER THE SAMPLE FREQUENCY” ,Samfreg 
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2300 GCLEAR 


231@ Tminc=!1./Samfreq 

2920 Timent=1Q@@@*#INT( Tmine/1.E-6) 
2339 TimeS="TIME "&VALS( Timent ) 
2342 Cat=VAL( Counts} 

293S2 Cokelvin=331.6 

2960 Num=VALCNS ) 

2979 Freqg=VAL( FreqS)*#1002 

2983 = 

ZOO) 

SCGO | Computation of the D-C offset prior to program run 
2e1id ! 

3320 =! 


3059 Off_ave=l 
53048 Of fset=2 
3052 FOR I=1 TO 19 


3360 CALL Read_ad200(Dat(*) Counts ,TimesS ,Off_ave) 
3079 FOR Ks) 00 (emt—)> 

3U39 Samave=SamavetDat(K ) 

53299 NEXT | 

35122 PRINT *€OuPUTING )-€ Geese. 


3112 NEXT I 

3129 Of f_ave=@ 

3139 Of fset=Samave/(1@*(Cnt-1)) 
3142 PRINT “OFFSET IS :” ,Offset 


Sisce 
Si6Ca 
las 1 computation of the speed of sound at a given temp and the 
3189 ! range of detection of the device 
3190 | 
3280 | 
2ci1@ Speed: INPUT “Enter the temperature (celsius) ° ,Temp 
3220 Spd=C@kelvin«(SQR( 1+ ( Temp /273))) 
3208 Rspd=Spd/z 
3240 Rdist=Tminc#Rspd 
5250 Maxrng=Cnt*Rdist 
3262 Exeng=((INT(Maxrng/32) )+.5)*230.@ 
sere Far=1./(Tminc*Spd) 
5282 Lambda=Spd/Fregq 
258 KZ=((2*PI)/Lambda)*(1/3) 
300 Pls_lng=(Rspd*Num)/Freq 
Sale CALL Attenuation(Freg Temp ,Atten?) 
os | 
Soae | 
3240 PRINT “Enter the relative intensity division levei. This 
3350 PRINT “value is used to determine the plot intensity by “ 
3362 PRINT “dividing the block average sum by this number. “ 
3578 PRINT “This value is dependent upon the gain of" 
53380 PRINT “the device and will probably need adjustment" 
YEG PRINT “during run. Start with a value of about 4000" 
3480 PRINT Pt 
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3412 
342 
243 


4 
344 


5450 
345 


aaa? 
ce re 


EvVimeuNPuT Ilva 
| 
| 
! set up the plot 
| 
{ 


Again:CALL Plot_setup(Nplot ,SitesS Maxrng ,Scale_y ,Exrng,Zone$) 
Nrec=0 
FOR I=! TO 202 
Ci( I= 
NEXT I 
Y=2 
Npoint1!=@ 
OUTPUT KBD; "L"3 
| 
| 
Sync: ! synchronize data collection with clock 
TI=INTCTIMEDATE MCD 864200) 
IF TI<TO THEN T@=TO-86400 
Pewitiewoc! (ten Gere Sync 
TQ=T I 
| 
t 
Data collection and reduction 
| 
| 
Read_sig: ! read the A-D converter 
CALL Read_ad2@@(Dat(*) ,CountS ,TimesS ,Off_ave) 
Nooint=INT(Ti MOD 368@ MOD Nplot) 
Npoint2=Npoint-Npointl 
Num] =Num- | 
Num2=Cnt-Num 
Ovd=!./(Num) 
K=Q 
| 
| 
FOR I=1 TO Nume STEP Num 
R=I*Rdist 
IF I<Far THEN R=] 
K=K + | 
Run_ave=@ 
FOR J=I TQ TiNunl 
X=Cait J )-Offset 
Run_ave=Run_aveti <A 
NEXT J 
Ns=SQR(Run_ave*Ovd ) 
Point(K )=R*(Ns-Corr ) 
IF Kl>1@ THEN 
CtCK aCtCK 40 (Point (K )%2 )*#EXPC2S*ReAtten) ) 
Comme CrCkh ets. o5se-c) 27 eer )/2imo 
ENO IF 
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3920 
3938 
3940 
Gaae 
3362 
2370 
3980 
3930 
4200 
4018 
4022 
4@30 
4040 
4050 
4062 
4372 
4080 
4090 
4100 
411@ 
4122 
4150 
414 
4152 
4160 
4170 
4180 
4198 
4200 
4212 


Bay 
a&& 


4230 
4242 
4252 
4260 
4270 
4280 
4290 
4390 
4512 
452 

4550 
43542 
4352 
4560 
4372 
4582 
4390 


IF Point(K )>ABS(32767) THEN Poimt(Ko~32767 
D2(K+1@ H=INTC Point (lt )) 
NEXT I 
IF K1>1@ THEN 
Y=Y+] 


Noise correction routine 


[TF Ns<Chimi tts) THEN 
Noise=NoisetNs 
Kl=Ki+1 
IF K1>1@ THEN 

Corr=Noise/Kl 
Limit=Corr 
ELSE 
Corr=0 
ENCe LE 
END IF 


| 
| 
| Plotting of the data 


Kmax=K 
REDIM O2(Kmaxt15) 
! positioning the data on the plot by time of trace 
| 
Time=(Npoint2/Nplot )*42@Q 
Timedist=2t( (Npoint2+Npoint! )/Nplot)*420 
IF Npointi=@ THEN Time=5 
IF Npoint<& THEN Timedist=6 
End_pt=Scale_y*#26@ 
D2¢(2 )=Kmax 
D2(3)=INT( Temp ) 
OY 4 =INT( Tine) 
D2¢(S)=INT(CTimedist+#1@) 
WINDOW @,428 2,260 
GRAPHICS ON 
{ 
! set the vertical increment 
Inc=( (Num#*Tminc*Rspd*tEnd_pt )/Maxrng) 
D2(6)=INT( Inc* 1800) 
| 
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Pp fF 
69 e3 0) 


fa 


compute the intensity of the return, move to the proper 
coordinates and plot the appropriate colored block 


| 
| 
| 
FOR K=! TO Kmax 
B=eoint (Kk )/livi 
IF Z>1 THEN Z=] 
IF Z<@ THEN Z=¢ 
AREA INTENSITY Z,Z,Z 
Dis=INT(K*Inc )+] 
MOVE Timedist ,Dis 
RECTANGLE Time,Inc ,FILL 
NEXT K 
| 
{ keep an account of the trace numbers taken on the plot 
| 
Nrec=Nrect! 
Rec_num=INT(Nrec? 
D2¢ 1 )=Rec_num 
Npoint!=Npoint 
| 


Save routine - if function key is set then 
the plot will be saved 


| 
! 
| 
| 
temcmecial) THEN 

ASSIGN @Filel TO FilelS 


OUTPUT @Filel ,Rec_num;D2(#) 
END IF 


| 
{ graphics dump of plot after 15 minute intervals 
{ 


IF Timedist>41S THEN 


Pet erp: PRUNE foc! 


Pein 7 

PRINT “ i 

PRINT “ ‘ 

peoemun=eLoemann | 

UMP GRAPHICS #781 

FOR I=! TO Kinax 
Gee =e ek CG Pema taio,) 2c) * G@iaeyecacd) thot Tey) 
Gol =6e (Loe C l/l 
Chior 1o8Cl e/ Ft) 
Ct(I)=Ct(1)/Pls_lng 
IF Ct(I)<1.@E-10@ THEN Ct( I )=1.0E-10@ 

NEXT I 

{ 


Sal 


| 

! Cl°2 Computations and pete 
| 

SCEEMI 

VIEWRORTei5, 120-15 ,89 

WINDOW @ 500 ,@ ,Exrng 
N\AGcwo0-Expnma/ 15.0 ,0 ,o0ee<ene, o 
CLIP OFF 

Col 4Abee oo 

LORG 6 

FOR M=0 TO 300 STEP 32 

MOVE Me -Eexrngs 45 

EABEL ~IeGE 3(<M/ 502410) 


NEXT M 

MOVE 14@,-18 
Coleewe 

LABEL LOG Oise rma 
LORG 8 


FOR M=@ TO sExgmo Sle INT Cexrmaec) 
MOVE O,M 
LABEL M 

NEXT M 

LOTR Try 2 

LORG 6G 

MOVE -40,Exrng/2 

Co pene 

LABEL "RANGE IN METERS" 

ED TRC 

LORG 4 

MOVE 150 ,Exrngt3 

LAB.) IMESAVERRSGED Cie "; Si ves 

CEleeo. 

MOVE 302 ,2 

FOR I=] TQ Kmax 
Ypl=(Maxrng/Kmax )#I 
Pnt=(LO6(Ct( 1))/2.302585 1) 
LeUP at 10 HEN ie no 
IF Pnt<@ THEN Pnt=2 
Pnt=Pnt *S@ 
DRAW Pat ,Ypl 

NE et 

DUMP GRAPHICS 4791 

PRINT < 


PRINTER IS CRT 
GCLEAR 

K1=12 
Limit=Corr 
Noise=!@*Corr 


Dic 


~J ~) -J 


orororeici c 
t) Co°=)°67 
C3 C2 C3 Cd OO 


-~j -J 


Ramen emi) Titi y 
ASSIGN @Filel TO * | 

END IF 

Cnirl=0 


e 


if the plot is to be saved, the file is created 
and named | 


! 
! 
! 
! 
IF Save_plt=1 THEN 
Covlmeverinit(Pise adcresss Nrec Filels) 
Save_plt=@ 
Cntrl=| 


END IF 

! 

! 

1! start the next 1S minute plot 
| 


GOTO Again 
Ener 


next trace 
toggle set if plot is to be saved 


| 
| 
| 
| 
GOTO Sync 


| 
Svpt: Save_plt=! 
GOTO Read_sig 


| 
|! completion routine 
( 
| 


-4 
=o 


Ou ris 


SUBROUTINE SECTION 


—-—- Se HMO 
= 
Oo 


cor 


6820 
6842 
62S 
6060 
6872 
6082 
6090 
6182 
6112 
G12 
6150 


61490 © 


6152 
6168 
6170 
6182 
61932 
6200 
' 6210 
6220 
G23 

6249 
6250 
6268 
fa 

6280 
6522 
6510 
6520 
Sage 
6542 
6552 
6560 
6579 
6580 
6530 
6602 
6610 
6620 
6630 


SUB Freq_init(FreqS NS) 
{ 


! 
f 
! 
| 


SETUP CF THE RP33i4A FUNCTION GENERATGR 


Ap$="AP" 
Fras="FR" 
NmsS="NM" 
En$="EN" 
Vos="V0" 
HzS="KZ"* 
INPUT “FREQUENCY DESIREO (kHZ) (S RECOMMENDED )” ,Freqs 
INPUT =AMPEITUDE BESIRED (V2) ssa amps 
IF VAL(Amp$)>1.S THEN 
Amps="{.5" 
PRINT “AMPLITUDE QF FUNCTION GENERATOR [5S 1.5 V" 
END IF 
INPUT “NUMBER OF CYCLES PER BURST CINTEGER) (18@ RECOMMENDED >)” ,NS 
CUTRUT 7o7--nes 
OUTPUT "7273 Ska 
OUTPUT 7O07;ApS&AmpS&VoSaFragS&FreqS&HzSaNmS8NS8ENS 


SUBEND 


| 
{ 
| 


SUB Init_adZz@@ 


| 
! INITIALIZATION GF THE A-O CONVERTER 
| 
] 


Ad_sel_code=!7 
Dummy=READIC. Ad_sel_code ,3) 
WRITCIO Ad_sel_ccode ,Q;2 
CONTROL Ad_sel_code ,Q; 1 


SUBEND 
| 
! 
| 
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7080 
7018 
7029 
7830 
7849 
7852 
7262 
7872 
7889 
7890 
7102 
ae 
7120 
noo 
7442 
7150 
7162 
7170 
7182 
7198 
20@ 
7210 
ee 
72350 
7240 
7258 
7260 
7279 


7288 SUBEND 


ado 
7500 
310 


SUB Set_time( Zones) 


. 
° 
* 
o 


! 
! 
! 


SET THE TIME DATE RECORDER 


eet ieewno i TIME REGERENGGeake YOU USING? INPUT: * 


PRINT “ ISGOR UNIVERSAL TIME 

Pa Ze rOR EOC ieee Line: 

PRINT “ 3 FOR YOUR OWN CLASSIFICATION" 
INPUT K 


IF K=2 THEN 
ZoneS="(LOCAL )" ° 
ELSE 
IF K#3 THEN 
INPUT “WHAT IS YOUR TIME REFERENCE" ,ZoneS 
Boe 
ZoneS="(UTC)” 
ENO IF 
ENO IF 
IF TIMEDATE<DATE( "14 AUG 1984") THEN 
INPUT “ENTER OB MMM YYYY° © Dates 
[INPUT “ENTER “em: MIN SSCs =  himes 
SET TIMEDATE DATE(DateS )+TIME( Times ) 
BRUNT CATES Th Pieenre (Ries ( PEMeBATE? 
Tstart=TIMEDATE 
TO=Tstart MOD 86400 
END IF 
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7509 
7510 
7520 
7530 
7542 
7550 
7560 
757@ 
7588 
75938 
7682 
7610 
7622 
7630 
7640 
7650 
7662 
767@ 
7682 
7693@ 
7782 
7710 
7720 
7730 
7742 
775@ 
7760 
7770 
778@ 
7732 
7800 
7310 
78a 

7859 
7940 
7350 
786@ 
7372 


SUB Read_ad2Q@@( INTEGER Dat (+) BUFFER ,CountS ,Times ,Of f_ave) 


| 
1 INFOTEK A-D ROUTINE SET UP FOR EXTERNAL TRIGGER 
| 
| 


Ad_sel_code=!17 
| 
{INITIALIZATION OF THE A-D CONVERTER 
QUTPUT Ad_sel_code; "RESET" 
QUTPUT Ad_sel_code;" INTERNAL” ,"“COUNT “&CountS, “HOLDON" 
QUTPUT Ad_sel_coda; "DELAYON“ “SELECT 1s! end“ ,Times 
OUTPUT Ad_sel_code; "STATUS" 
ENTER Ad_sel_code;Resp$ 
IF Resp$=*-------- “ THEN 
ASSIGN @Ad20@ TO Ad_sel_code;WORD 
| 
| triggering of the HP3341A 
IF Off_ave=0 THEN 
TRIGGER 707 
ENO IF 
ASSIGN @Buf TO BUFFER Dat(*) 
TRANSFER @Ad2@@ TO @Buf; WAIT 
QUTPUT Ad_sel_code;" “ 
OUTPUT Ad_sel_code; "STATUS" 
ENTER Ad_sel_code;RespS 


IF Resas<7  -------= " THEN — 
PRINT "“ERROR= “;Resp$ 
END IF 
BLS = ; 
PRINT “ERROR DURING INITIALIZATION =";Resps 
ENiaeis 
SUGEND 


! 
! 
! 
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8822 
8212 
8020 
8832 
8242 
8052 
8962 
3879 
8880 
8092 
8182 
8112 
8122 
8130 
8142 
8152 
8162 
8172 
8180 
8190 
8200 
8218 
ee 
8230 
8242 
8250 
826@ 
8272 
8230 
8292 
8300 
8312 
8320 
8322 
8340 
8552 
8362 
837 
8595 
8530 


ATA 
wud 


8410 
8422 
8450 
8442 
8450 
8462 
8472 
8432 


SUB Plot_setup(Nplot ,SiteS Maxrng ,Scale_y ,Exrng ,ZoneS) 


oes OF The Tine Peol ON THE CRI 


Scale_y=Maxrng/Exrng 
GRAPHICS ON 

LINE TYPE 1} 

VIEWPORT 15,120,15,8@ 
WINDOW @ ,Nplot ,@,Exrng 
AWAGoeNplot? Wal Exrng/|G.0-0 Nplot/3 ,exrng/s 
CLIP OFF 

CSIZE 4,.6 

LORG 6 

TI=TIMEDATE MOD 8642@@ 
Hres=T! DIV 36982 

T2=T! MOO 3600 

Min=T2 DIV 62 


Qtrhr=Min DIV 1S 


FOR M=@ TO Noelot STEP INT(Nplot/3) 
MOVE M.—exrng/45 
Qtrmin=Qtrhr#ttStCMeZ/Nolot #5 
IF Qtrmin=6@ THEN 

Qtrmin=9 
Hrs=Hrstl 
END IF 
LABEL USING °00,A8,ZZ°3Hrs;3°:°;Qtrmin 

NEXT M 

Meve Nelot/2 ,-15 

EABEL “TIME “&Zones 


LABEL ORDINATE 


LORG 8 

FOR M=@ TO Exrng STEP INT(CEx«rng/3) 
MOVE @,M 
paises im! 

NEXT ™ 

Eine rife 

LORG 6 

MOUE =Nealot/?7.Exrno/ed 

LABEL "RANGE IN METERS" 


ions 


LOIR @ 
LORG 4 
MOVE Nplot/2 ,Exrngt3 


a 


8490 
8500 
851@ 
8520 
8550 
8540 
9209 
9310 
9020 
9052 
9242 
90S@ 
J@62 
3070 
3@82 
9090 
3182 
9112 
9122 
9130 
9142 
9158 
3162 
317 
3180 
919@ 
3280 
9210 


aed 
ha he 


3250 
9249 

292 
3260 
9272 
9280 
J29@ 
9582 
931@ 


“=. {2% 


9549 
J250 
3369 
937@ 
3380 
3390 
9422 
9412 


LABEL DATES TIMEDATE a = ";Sites: 


CLIP ON 
SUBEND 
UB Attenuation( Freq ,Temp ,Atten) 
This subprogram calculates the attenuation of the 
sound in air based upon equations in Neff 1975 
(source of subroutine: Thesis of R. Fuller) 
Variables 
Atom_pres - input atmospheric pressure in mb 
Atten - attenuation coefficient of acoustic wave 


at the frequency of the maximum attenuation 
for the given input conditions. 


F - Ratio of the frequency to frequency at maximum 
attenuation. 


Fmax —- Frequency of the maximum attenuation 

H - variable used in the integration of excess attenuation 
Pstar - variable used in intermediate calculations 

Tstar - variable used as an intermediate in calculation 


of attenuation. 


| 
| 

S 

! 

| 

| 

| 

! 

! Att_max - Variable in program. It is the attenuation 
| 

| 

! 

| 

| 

| 

| 

| 

| 

| 

! Wat_pres - Atmospheric water pressure in mb. 
INPUT “Enter the atmospheric pressure in mb" ,Atom_pres 
INPUT “Enter the atmospheric water pressure in mb" ,Wat_pres 
H=1OG*Wat_pres/Atom_pres 

Tstar=(1.8*Tempt+492 0/519 

Pstar=Atom_pres/1@1i4 

Fmax=( 1@+6600*H+44400*#H*#H )*Pstar/Tstar”*.3 
Att_max=.@078*Fmax*Tstar*(-Z.5 )*EXP(7.77#(1-1/Tstar)) 
F=Freq/Fmax 
Atten=(Att_max/304.3)#¢( (. (8"F }°2+(2aFeF/( THF *#F 7297.5 
Atten=(Attent!.74E-1@#FreqtFreg )/4.35 
SUBEND 
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9420 
9430 
9582 
9510 
ace 
95590 
9542 
JSS 
9560 
9572 
9580 
9590 
9602 
361@ 
9620 


SUB File_init(Disc_address$ ,Nrec ,FilelS) 


! 
 Semente Views lORAGCE FILE ON THE DIS€ FOR 
| THE REDUCED DATA 

| 

| 


Theo) en leER VHe REQUCEO DATA OUTPUT FILENAME “Files 
FilelS=FileS&Disc_addresss 
CREATE BDAT File!$,180,400 ' 
ASSIGN @Filel TO Filels 
SUBEND 


og 


APPENDIX B 


SPEAKER AND ARRAY ANALYSIS 

In the design of our echosounder, it was determined 
that a rapid decay time was required to obtain accurate 
short range information. Weight restrictions involved with 
equipment transportation plus maximum response at high 
frequencies led to our decision to use piezo ceramic 
speakers. The Motorola KSN 1005A speaker was selected 
based on these requirements and the specifications charted 
in the Motorola catalog [Ref. 16] and reproduced below in 
Figures 23 through 25. 


Nominal Power, Impedance, and Distortion Ratings 
a a a 





MAXIMUM ALLOWABLE 
SPEAKER (RMS) 


VOLTAGE ACROSS 


FREQUENCY (Hz) 
300 500 1000 892000 5000 10000 20000 40000 


Fig. 23. Speaker Ratings 
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Typical Frequency Response 
Loki] Toi FF igh hgh UG hUTGlhUGlLUT UT] FT UF OUT )hUYlhlUWPC CUM 
Ley | fcr | ese | sees | es | | | a en ge | 
= 
INPUT. 2.8V 
MICROPHONE DISTANCE: 1/2 Meter 
(1 Watt tnto 8 Ohms) 


OUTPUT (dB) 





FREQUENCY (Hz) 
500 1000 2000 §000 10000 20000 40000 


Fig. 24. Speaker Frequency Response 





Dimensions: KSN 1005A, KSN 1003A 


ER Oy TOT ME ow te 


FOUR 218" 15 Smm) DIAMETER HOLES EQUALLY 
SPACED ON A 3.947 (100,1mm) OIAMETER B.C. 





d 
j y 

¥ 334 

y = 

——— | 

269 / ( ie ee ie 

Pad \/ | 
070. \ Q pL 


WEIGHT: 75 GRAMS oe a us , ; 


Fig. 25. Speaker Dimensions 


oe 


The frequency response charts indicate that a maximum 
response for our speakers occurs at a resonant frequency of 
5000 Hertz. This frequency was used as the baseline from 
which all our measurements are made. 

Using the speakers in an anechoic chamber, the average 
e-1 voltage decay time was measured to be approximately 900 
usec (Fig. 26). This decay time Peaneiates into a sound 
propagation distance of just over 15.0 centimeters (at STP) 
from the speakers. Considering our requirements, this 


speaker is ideally suited to serve our purpose. 


At=2,.924 ms 
At Av= 25 mv 


Wha. NANA ARAA 








Fig. 26. Speaker Decay Time Trace 
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pe ACOUSTIC ARRAY 

Our next consideration was the echosounder array 
pattern. Ideally, the acoustic sources should be placed 
exactly one half wavelength apart. At a frequency of 5000 
Hertz, this would require spacings of 3.4 centimeters (at 
STP) which is physically impossible for the speakers we 
have chosen. The closest possible Spacing is 7.62 
centimeters between sources after shaving off the flange of 
the horn (Fig. 25). 


From Kinsler, et al. fRef. 14], the equation for the 





directionality factor of a simple line array is derived as: 


N 
Sin { --- kd sin 6 ) 
Ji 2 
HCN 6) = eee ee ee . 
N i 
Sin ( --- kd sin @ ) 
Z 
where 


k is the wave number (27/ XX ), 

dis the distance between sources, 

N is the number of sources, and 

8 is the angle measured from a line perpendicular to 
the array to the direction of interest. 
However, this equation assumes simple point sources which 
does not adequately describe our speakers. It was necessary 


to couple this equation to the directional factor fora 
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piston source which is also identified in Kinsler, et al. 


[Ref. 14: p. 108] as: 


D(@) = | -------------------------- . 


where 

k is still the wave number, 

a is the radius of the piston source, 

8 is the angle measured from a line perpendicular to 
the array to the direction of interest, and 

Ji ( x ) ais a Bessel function of Order 1 with an 
argument x. 

These two equations were combined together to form the 
equation of directionality for a linear array of piston 
sources, L(N,8), by simple multiplication. 

iCNeie ee = Ne Oe eee) 

This equation was incorporated into a computer program, 
which may be found in Reference 3, and used to generate 
linear array patterns for various numbers of piston sources 
(Figs. 27 throucn 3b)2 

It was concluded from these plots that five linear 
elements would provide the best combination of forward 
directionality, sidelobe suppression, physical size anda 
relatively low cost. Then, in order to enhance both array 
efficiency and symmetry, we settled on a five by five 


element array design with 7.26 centimeter spacing between 
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PATTERN FOR 3 ELEMENTS PATTERN FOR 4 ELEMENTS 


Q 8 
° oO 
= a 
im 20 a 
rf O 
3 z 








J 
ud 
q 


Fig. 27. Three Element Array Fig. 28. Four Element Array 


PRTTERN FOR 5S ELEMENTS 


= 
© 






DECIBELS 
nu 
q 


Fig. 29. Five Element Array 


PATTERN FOR 6S ELEMENTS 9 PATTERN FOR 27 ELEMENTS 


-19 -12 






Wn 
g d 
20 bas 
-39 YW fr 
Wy 
Fig. 30. Six Element Array Fig. 31. Seven Element Array 
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speakers in both the vertical and horizontal directions. 
After verifying the manufacturers polarity designation 
for 35 speakers, we obtained Lissajous plots for each 
individual speaker (Fig. 32). Based on the speaker’s output 
to input voltage ratios as illustrated by these plots, we 
were able to rank all our speakers by signal efficiency. 
It was based upon this criteria that we selected the 25 
most efficient speakers for the array, placing the best 
speakers at the center and subsequent ranking speakers 


further toward the sides and corners. 






SS say 
spea 

steate, ee | er 
Pace tT 
|| btretertnce sheaker|inpurt——ol| 



















Fig. 32. Lissajous Plot 
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The 25 selected speakers were mounted in a five by five 
planar array on a balsa wood insulated bilayered sheet 
metal board. After wiring all the speakers in parallel, we 
surrounded all the electrical connections and speaker backs 
with two 3.0 centimeter layers of foam insulation 
sandwiching a 1.0 millimeter lead sheet. Then the entire 
array mounting was enclosed in a 44 by 44 by 5 centimeter 
sheet metal box. This design was chosen to suppress 
virtually all acoustic energy propagating out the rear 
hemisphere of the array, while shielding the array from any 


external electrical interference (Figs. 33 and 34). 
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Fig. 33. Array Photo 
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Array Photo 


34. 


Fig. 
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APPENDIX C 


ENCLOSURE DESIGN 


Acoustic echosounding has proven to be an extremely 
useful technique for probing and analyzing the lower 
atmosphere. In order to most efficiently utilize the 
acoustic waves transmitted and later received by this 
remote sensing method, it is essential to have an efficient 
antenna with highly directive beams and strongly suppressed 
sidelobes. Antenna design becomes increasingly more 
important in a noisy environment where noise pollution 
Within the sidelobes may dominate the desired signal within 
the main lobe. Hall and Wescott [{Ref. 12] showed that 
Sidelobe suppression improved with higher frequencies. 
Their studies showed that the measured 90 degree sidelobe 
Suppression ranged from 38 dB at 1 KHz to 50 dB at 5 KHa2. 
Furthermore, any significant improvement in sidelobe 
suppression could only be obtained by surrounding the 
antenna with an acoustic energy absorbing cuff or shroud. 

In an effort to maximize our antenna main lobe to 
Sidelobe power ratio, we intend to operate only at high 
frequencies as discussed in the previous section. 
Additionally, we have designed an acoustic energy absorbing 


enclosure. Many designs were considered based upon 


eo 


previous research in the field of echosounding [Refs. 17 
through 20). In addition, we obtained the actual acoustic 
beam patterns for our array uSing a computer program 
written by LCDR Butler [Ref. 21] which we eae T Orgone 
purposes. This modified version of LCDR Butler’s program 
may be found in Reference 3. By rotating the array in an 
anechoic chamber, we were able to produce highly accurate 


polar plots of the array beam patterns (Figs. 35 and 36). 


5 X 5 ACOUSTIC ARRAY BEAM PATTERN 





R=3.60 M MIKE VT= 8.2465 INPUT=5.0 V FREQ=5. 0 KHZ 


Fig. 35. Polar Plot of 5 X 5 Array Beam Pattern 
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DIAGONAL 5 X 3 ACOUSTIC ARRAY BEAM PATTERN 
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Fig. 36. Polar Plot of 5 X 5 Array Beam Pattern 
at a 45° Aspect 
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These polar plots and the actual coordinates 
corresponding to the individual data points confirmed the 
computer prediction for the five element line array (Fig. 
29). Additional polar plots obtained by varying both the 
input array voltage and range between array and microphone 
further support our claim that for a 5 KHz carrier 
frequency, the main lobe is confined to a divergence angle 
of 20 degrees. since it is our aim to suppress all 
sidelobes and utilize solely the main lobe, we chose not to 
taper our enclosure as most previous researchers had. 
Rather we designed the enclosure based upon the dimensions 
of the array itself and the acoustic beams it generated. 

Plywood was used for the construction of the enclosure 


and provided not only a rigid, inexpensive framework, but 


also proved ce greatly attenuate external noise 
interference. Anticipating all kinds of weather conditions 
during data collection, the plywood enclosure was first 


waterproofed with four coats of marine varnish. Grooved 
joints, caulking and weather stripping were also design 
considerations. 

A millimeter layer of lead can suppress an acoustic 
Signal as much as 40 dB (Figs. 37 and 38). About 7.0 
centimeters of corrugated, egg-carton design foam can 
suppress a signal another 3 to 4 dB (Figs. 37 and 39). 
Together they make an extremely efficient absorbing 


material for use in our enclosure. 
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Fig. 37. Response Reference, No Insulation 
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Fig. 38. Signal Suppression by Lead 
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Fig. 39. Signal Suppression by Foam 


Lead-lined absorbing foam was glued to all inner 


surfaces of the enclosure with two 1.0 millimeter layers of 


lead overlapping at all corners. Strong aluminum brackets 


were used to connect the four side panels to each other as 


well as to the enclosure base (Figs. 40 through 42). 
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Fig. 40. Enclosure Photo, Fully Assembled 
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Interior 


Enclosure Photo, 
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ECHOSOUNDER OUTPUT 

Sixteen echosounder output traces were included to 
exhibit typical atmospheric activity. Mam, Of the plots, 
such as the 14:45, 15:30, and the 17:30 of 26 April and the 
TO: OGRE. to eee OC 4tebe Toei bemeeooi5, and the 17:15 of 
2? April, have convective plumes which are prevalent 
whenever a heat flux between the surface and atmosphere 
exists. Hhew@plOovuSs Of !Os1o ana £9:15 on the 26th of April 
and of 18:00, 18:45 and 20:00 on the 27th of April clearly 
show the passage of the neutral event which is encountered 
when the atmospheric and surface temperature difference 
becomes negligible. Finally, the plot of 16:30 on 26 April 
can be associated with strong winds which exhibit a 
somewhat uniform return for all altitudes across the entire 
15 minute interval of the trace. 

The horizontal axis labels of the Cr2 plots are small 
and may be hard to. read. This axis is a log scale 
beginning with 1°10-1°9 at the left side of the plot and 
ending with 1°10° at the right side of the plot. Each tick 
mark moving left to right along this axis represents an 


integer increase in the exponent of 10. 
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